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ABSTRACT 

^' 

Q ! Detection of ^Li has been shown for energetic solar events, one chromospher- 

^ ! ically active binary, and several dwarf halo stars. We had previously found a 77^ 

ci , = 0.03±0.01 for active K dwarf GJ 117 using VLT UVES observations. Here we 



> 



present high signal-to-noise (>1000) high spectral resolution observations taken 
with the McDonald Observatory's 2.7ni and echelle spectrometer of GJ 117. We 



cn ! have used the solar spectrum and template stars to eliminate possible blends, such 

\^ [ as Ti I, in the ^Li spectral region. Our new analysis, using an updated PHOENIX 

^ I model atmosphere finds 7^ = 0.05±0.02. Additionally, bisector analysis showed 

^ I no significant red asymmetries that would affect the lithium line profile. No 

CX) ! changes above the statistical uncertainties are found between the VLT and Mc- 

^-^ ' Donald data. The amount of ^Li derived for GJ 117 is consistent with creation 
in spallation reactions on the stellar surface, but we caution that uncertainties 

rN I in the continuum level may cause additional uncertainty in the ^Li fraction. 

^'. 
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Introduction 



The destruction of primordial lithium in low-mass stars occurs during the pre-main- 
sequence state and little or no Li is therefore expected in these objects when they arrive on 
the main sequence. Studies of open clusters found th at the largest lithium abundances occur 
for the fastest rotating stars (jSoderblom et al.lll993l ). This result is in contradiction to the 
expectation that fast rotation enhances the mixing process that leads to increased lithium 
depletion, and the rotational history must play an important role in the mixing process 
( jPinsonneaultl 1 1 99 71 ) . 



One of the key assumptions in the models of Li depletion, is that lithium can not be 
produced on the surface of late-type stars and its abundance will only decrease as a function 
of time. However, if lithium can be produced on the stellar surface, a " quasi-equilibrium" 
abundance would be reached which could depend on the activity level. Thus, we have been 
investigating if ^Li can be formed in active stars by spallation in stellar flares as a result of 
their high levels of activity. 

The presence of ^Li in metal poor halo stars sparked much debate in the early 1990s 
( ISmith et al.lll993l ). Subsequent observations confirmed the detection of ^Li in such halo 
stars and its production was attributed to f ormation by cosmic ray spallation in the ISM 
( Smith et al.lll998l: Hobbs &: ThorburnI 119971 ). although more recent work argues for a pre- 



Galactic origin (JAsplund et al 



200 
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The detection of ^Li en hancement has been repo rted during an energetic flare on a chro- 
mosp herically active binary (JMontes fc Ramseylll998l ). and an active K dwarf (IChristian et al. 
2005 ). There are also indicat ions that this isotope may be enhanced in sunspots (JGiampapa 
1984J : iRitzenhoff et al.l 119971 ). The ge neration of significant ^Li in stellar and solar flare s 
has been predicted by several authors (ICanal et al.lll975l : IWalker et al.lll985l : iLivshitd 119971 ). 
Such enhanced ^Li a bundances should also be ex pected for stellar objects if extreme energetic 
conditions are met (jsee MuUan fc Linskylll999l ). 



In an earlier paper we have shown the de tection of ^Li in the a tmosphere of the active 
K dwarf GJ 117 (=HD 17925) at the 3% level Jchristian et al.ll2005h . Given the signiflcance 
of this result we have decided to re-observe GJ 117 at a much higher signal-to-noise than 
observed so far . Our observations, analysis, and models (which are given in more detail in 
Christian et al.l ( 20051 )) are presented in § 2. Our results are discussed in § 3 with concluding 
remarks presented in § 4. 
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2. Observations & Data Analysis 



2.1. Observations 



Observations of GJ 117 were conducted between 03 - 09 November 2006 using 107 inch 
Harlan J. Smith telescope and echelle spectrograph at McDonald Observatory. The CS21- 
el instrument was used with the El grating centered at 6708A and the TK3 2kx2k CCD. 
This set-up provided a resolution of ~130,000. Multiple exposure of GJ 117 were taken for 
a combined spectrum with 55 ksec of exposure and a signal-to-noise of over «il300. The 
data were reduced with standard IRAF tasks, imred.ccdred and twodspec.apextract, 
and further analysis carried out using IDL and the starlink based DIPSO software. The 
spectrum of GJ 117 compared to the solar spectrum (taken with the same set-up) in the Li I 
6708A region is shown in Figure [H 

We also obtained Extreme Ultraviolet Explorer (EUVE) data from the Multimission 
Archive at Space Telescope. EUVE observed GJ 117 for nearly 200 ksec in December 1994 
with a mean count rate in the Deep Survey (DS) 100 A bandpass of 0.09 counts s~^. We 
reduced the EUVE data with the standard IRAF EUV reduction tasks and constructed 
light curves for the D S data using the xray.xtiming package as discussed previously in 
Christian et al.l (119981 ). A log of the GJ 117 observations are listed in Table 1. 



2.2. Model Spectra 

^Li is separated from the ^Li doublet only by 0.16 A and has a strength typically less 
than a few percent of that feature. In addition to the feature's intrinsic weakness, the 
situation is further complicated by b lending with nearby lines, such as CN, Fe I, and Ti I 
( iNissen et al.lll999l : iReddy et al.ll2002l ). Contamination by such lines becomes more prevalent 
as we move down the main sequence to cooler metal-rich stars. Template stars were observed 
along with GJ 117 to try to resolve the issue of blending. These stars were chosen to have 
a log(T) similar to GJ 117, and of the several observed, GJ 211 ( HD 37394) also has t he 
closest metallicity (Fe/H = 0.12), compared to GJ 117 (;Fe/H=0.15: lluck fc Heite] J2005h l 
The spectra of GJ 211 and GJ 117 are compared in Figure [2l CN lines near 6707.5 that 
have been seen to be < 1 niA in the solar spectrum, are also very weak in GJ 211, with an 
equivalent width upper limit of ~0.5 mA. CN should have a negligible effect on the GJ 117 
Li I line profile. Additionally there is no evidence for a strong Ti I line near 6708.025 A with 
a conservative upper limit of < 1 mA. 



The general stellar atmosphere code Phoenix ( see Hauschildt et al.lll995l : lAllard fc Hauschildt 
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Table 1. OBSERVATION LOG 



Name 


Observation 


Exposure 


Air Mass 




Date 


(sec) 






McDonald 




GJ 117^ 


2006 Nov 04 


1800 


1.44 






1800 


1.39 






1800 


1.38 






1800 


1.39 






1800 


1.56 




2006 Nov 05 


1800 


1.45 






2500 


1.39 






2500 


1.38 






2000 


1.65 




2006 Nov 06 


1800 


1.49 






3000 


1.42 






3000 


1.38 




2006 Nov 07 


3000 


1.39 






1800 


1.48 




2006 Nov 08 


1800 


1.68 






3000 


1.52 






1800 


1.86 




2006 Nov 09 


1800 


1.42 






3000 


1.38 






1800 


1.46 


GJ211'^ 


2006 Nov 04 


900 


1.13 






1800 


1.65 




2006 Nov 05 


900 


1.10 






900 


1.11 




2006 Nov 06 


1200 


1.15 




EUVE 




GJ 117 


1994 Dec 02 


192^ 





'^GJ 117: RAj2ooo= 02 52 32.1; Decj2ooo= -12 
46 11.0 

''Stellar template (see text). 

'^EUVE observation time in kiloseconds. 
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Fig. 1. — McDonald echelle spectra of GJ 117 (solid histogram) compared to the solar 
spectrum (dashed line) in the Li I 6707. 8A region. The spectra are shown in total counts 
and the solar spectrum is offset by 30% for ease of presentation. 
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Fig. 2. — Comparison of the McDonald GJ 117 spectrum (solid histogram) to that of a K2 
template star, GJ 211 (dashed line). 
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19951 ) was used to calculate the oretical stellar spectra. We have us ed the Phoenix NextGen 
series LTE atmosphere models (JHauschildt. AUard. &: Baronlll999f) for the effective tempera - 
ture (log(Te//)=3.7) and gravitv flogfq)=4.6) of GJ 117 JCravlbood : Icincunegui et al.l[2007l ). 
We constructed a grid of models fro m NextGen that sampl ed every 200K in Teff and O.Sdex 
in \og{g). Direct opacity sampling (JHauschildt et al.ll200ll ) is used within Phoenix to han- 
dle any line blanketing, where the opacity at each wavelength point is calculated as a sum 
of opacities from all the contributing s pecies. We added the ^Li resonance do ublet to the 



Phoe nix master line list (JKuruczlll995l ) using the wavelengt hs and qf valu e s of 



Smith et al. 



( 19931 ). A Ti abundance of log(Ti) = 5.22 was chosen from lLuck fc Heiterl ( 20051 ). We 



cau- 



tion that the Ti abundance and jj^ ratio are anti-correlated and increasing the Ti abundance 
decreases the jj^ ratio and conversely, decreasing the Ti abundance increases the y^ ratio. 

Model spectra were calculated with a microturbulence velocity of ^=1.5 km s~^. This 
value is ty pical for solar-type s t ars a nd is consistent with the microturbulence given for 
GJ 117 in kuende Prieto et al.l J2004h . We derived a ^sir ii of GJ 117 of 5.5± 0.2 km s"^ 
using x^ analysis of several lines in the Li I 6707. 8A region ( Christian et al.ll2005l ). We show 



model fits for the Fe I A6713.7 and Ca I at A6717.685 line profiles in Figure [3] using the 
derived fsini of 5.5 km s~^. 



Reddy et al.l (|2002| ) have shown that the apparent detection of ^Li for planet hosting 
star HD 82943 could be explained by the presence of Ti I lines near 6708.03 and 6708.1 A. 
Although our template spectrum of GJ 211 rules out any strong Ti I features in the red-side 
of the lithiur a line profile we cori structed a second set of Phoenix models (PHX2) with the 
Ti I values of iReddy et al.l (120021 ) for comparison and completeness. We show t he line list for 



i mpor tant lines in the vicinity of the lithium line profile and their comparison to lReddy et al. 



( I2OO2I ) in Table 2. This second Phoenix model (PH. X2) has the s a me va lues as shown for 
the first Phoenix model in Table 2, but includes the iReddy et al.l (120021 ) values for the Ti 
lines at 6707.752, 6708.025 and 6708.125 A. 

Bisector analysis has been shown to be useful in disentangling subtle enhancements in 
line profiles that may be caused by blending with weaker line s and indicate othe r effects in 
the stellar atmosphere, such as granulation. In a recent paper, ICayrel et al.l (120071 ) have used 
bisector analysis in observed line profiles of HD 74000 and have shown that previously un- 
detected line asymmetries may cause derived lithium fractions to be lower than previously 
reported. We have performed bisector analysis on 3 Fe lines (6703.5, 6705.1, 6713. 7A) in 
the Li 6708A region to evaluate the effect of asymmetries. We show the line bisectors for 2 
of these Fe lines in Figure JH We find that the line bisectors for these lines on GJ 117 are 
vertical, and show a blue asymmetry near the top of the profile, but show no red-symmetry. 
We estimate an upper limit to red-asymmetry of Ri0.005A (0.2 km/s). This is consistent 
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Table 2. Lines in the vicinity of the Li I 6707.8 A profile 



Wavelength Element LEP (eV) 



log fff 



A 




^2=^ 


PHX^ 


R02 


PHX 


6707.3810 


CN 


1.83 


1.83 


-2.170 


-2.141 


6707.4330 


Fel 


4.61 


4.61 


-2.283 


-2.283 


6707.4500 


Smll 


0.93 


0.93 


-1.040 


-2.379 


6707.4640 


CN 


0.79 


0.78 


-3.012 


-2.921 


6707.5210 


CN 


2.17 


2.00 


-1.428 


-1.358 


6707.5290 


CN 


0.96 


0.96 


-1.609 


-1.594 


6707.5630 


VI 


2.74 


2.74 


-1.530 


-1.995 


6707.6440 


Cr I 


4.21 


4.21 


-2.140 


-2.667 


6707.7520 


Til 


4.05 


4.77 


-2.654 


-2.656 


6707.7610^ 


^Li 


0.00 


0.00 


-0.002 


-0.009 


6707.7710 


Cal 


5.80 


5.79 


-4.015 


-4.014 


6707.8160 


CN 


1.21 


1.21 


-2.317 


-1.967 


6707.9130^ 


^Li 


0.00 


0.00 


-0.807 


-0.309 


6707.9210 


^Li 


0.00 


0.00 


-0.002 


-0.005 


6708.0250 


Til 


1.88 




-2.252 




6708.0728 


6Li 


0.00 


0.00 


-0.303 


-0.309 


6708.0940 


VI 


1.22 


1.22 


-3.113 


-3.113 


6708.1250 


Til 


1.88 


5.65 


-2.886 


-2.801 


6708.3750 


CN 


2.10 


1.89 


-2.252 


-1.959 



^R02 - from the line list of iReddv et all (l2002f ): PHX - from 
the PHOENIX line list used in this work. 
^Ave of AA6707.754 and 6707.766 ^Li lines 

=Ave of AA6707.904 and 6707.917 ^Li lines 



D 



0.9 - 




0.8 



6713.4 6713.7 



1.0 


Ca 


1 , , , 1 


0.8 




/ 


0.6 






0. 




V 

1 1 1 1 1 



5714.0 6717.2 6717.6 6718.0 

o 

Wavelength (A) 



Fig. 3. — Sample synthetic line profiles for: Fe I A6713.7 (left), and Ca I A6717.68 (right) 
lines using the derived vsini of 5.5 km s~^ (see text). 
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with our earlier estimates from VLT UVES spectra of GJ 117. We therefore conclude, that 
there are no significant red asymmetries that may affect the Li line profile. We discuss 
uncertainties in continuum placement in Section 3.1, which compares the McDonald and 
VLT spectra of GJ 117. 



3. Results & Discussion 

3.1. The 7^4 isotope ratio 

We have compared the combined McDonald spectrum to the Phoenix models using 
a least square fitting technique where t he quality of the fit is determined by x^ statistics 



( IGhristian et all 1200.4 ISmith et al.lll998f ). 

^' dof^ af ^'^ 

where Di and Mj are the data and model fluxes at data point i, respectively. We defined 
the variance, ctj as the square root of the counts at each wavelength, a ranged from ^ 900 in 
the lithium line center to Ril350 in the continuum. The standard continuum normalization 
was determined using the DIPSO CREGS routine for continuum placement. We performed a 
least-squares fitting to check this continuum level. We selected 0.3 A regions of continuum to 
the blue and red side the ^Li line profile. We fitted the model against the data and adjusting 
the data by a normalization factor between 0.98 and 1.02 in 0.001 increments, and determined 
a minimum x^ for a normalization factor of 1.003. We then used this normalization factor in 
fitting each PHOENIX ^Li model to the data by varying the 7^ fraction in 0.01 increments. 
In this way the most probable model with x^ ~1 was determined, x^ was computed over the 
red-side of the ^Li + ^Li line profile in the wavelength range of 6707.65 to 6708. 35A, giving 
54 degrees of freedom (dof). We show the combined McDonald spectrum and its comparison 
to 3 models with the 7^ fraction ranging from 0.0, 0.05, and 0.10 in Figure [5^, and an 
expanded view of the red-side of the lithium line profile in Figure Eb- The best fitting model 
had a 7^ = 0.05 with a x^ oi ~58 and the model with no ^Li had a x^ ^ 284. We show 
the reduced x^ values for each model and Table 3. For the residuals Figures [S^b we include 
a sign factor to x^ to indicate whether the observational data was larger than (positive) or 
smaller than (negative) the model. 

We computed the F-statistic and compared the reduced x^ for each model to the model 
with the minimum reduced x^ (model with 7^ = 0.05), and show this in Figure [6l For 54 
degree of freedom the 90% confidence limit for the F ratio is ~ 2.8 and this ratio is also 
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Fig. 4. — Sample line bisectors for 2 Fe I lines (A6703.5 & A6705.1) in the vicinity of the 
lithium line profile. 
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shown in Figure [6l The la uncertainty of the ®Li fraction is 0.02, giving -ff-. = 0.05±0.02. 



Recently, ICayrel et aLl (120071 ) cautioned that the continuum placement and the uncer- 
tainty in the continuum level may dominate the errors in the ^Li fraction, and we tested 
the effect on the this level on the ^Li fraction. We find that a normalization of 1.01 (1% 
increase in the continuum level derived a 7^ = 0.03 when re-fitting the lithium line, but 
clearly does not fit the line center (data now higher than the model) and this fit has x^ 
about 60 higher than the best fit normalization of the continuum, which gives 7^ = 0.05. 
Conversely, by lowering the continuum fraction by 1%, we find a txi = 0.06, but again this 
shows an overall poor fit to the line and continuum and has x^ about 100 higher than the 
best fit normalization of the continuum. 

We also investigated if small wavelength shifts could change the derived 7^ fraction. 
Our PHOENIX models are generated with a wavelength increment of 0.01 A and then re- 
binned and interpolated onto the observed wavelength scale. Thus any small shifts caused 
by an uncertainty in the wavelength scale should be accounted for as the data and model 
shift relative to each other. However, we tested how small absolute shifts between the data 
and model would change x^ aiid the derived y^. Shifts to the data relative to the model 
on the order of 1% of a resoluton element (0.5 mA) only changed x^ by a small amount 
and caused no changes to the computed jjj. However, we then tested what amount of an 
absolute shift would cause a change of 0.01 in the 7^. A shift of +5mA (10% of a resolution 
element) is needed to shift the best fitting model for the 7^ to 0.06 and similarly a shift 
of -5 mA gives the best fitting model to 0.04. However, both of these shifts produce a x^ 
approximately 30 higher than the best fitting model, and such a large shift is not supported 
by the wavelength calibration and modelling, and any smaller shifts are accounted for in the 
current uncertainty. 

We have found a 7^4 of 0.05±0.02 for our McDonald observations of GJ 117. This 
value is ~3% higher than our earlier VLT estimate. The main difference can be attributed 
to not i ncluding the s trong Ti I line at 6708.025 A. If we include this line at the strength 



used by iReddv et al.l fl2002h we find a 7^ of 0.02±0.01. The fitting results for this model 
are also summarized in Table 3. Our template star, GJ 211 and other K2 stars in our 
sample do not show any additional absorption near this wavelength stronger than 1 mA and 
our original PHOENIX model appears to be more consistent with GJ 117. An analysis of 
the VLT spectrum with our first PHOENIX model also finds a similar result for the 7^ of 
0.05±0.03 and is thus consistent with our higher signal-to-noise McDonald data. The spectra 
from the two epochs do not show a significant change above the statistical uncertainties and 
uncertainty in normalizing the spectra. We note that the VLT spectra have a S/N of ~ 400, 
significantly lower than the McDonald observations. 
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Table 3. Results for Xr for PHOENIX model fitting 



6 Li 


PHX'^ 


X'r 

PHX2^ 


0.00 


5.27 


1.70 


0.01 


3.52 


1.18 


0.02 


2.29 


0.98 


0.03 


1.56 


1.11 


0.04 


1.14 


1.45 


0.05 


1.08 


2.06 


0.06 


1.42 


3.01 


0.07 


2.15 


4.01 


0.08 


3.01 


5.36 


0.09 


4.42 


6.88 


0.10 


6.09 


8.76 



'"PHX - from the 
PHOENIX line list used in 
this work. 

^PHX2 model with addi- 
tional Ti I from the line list 
of Reddy et al. 2002 
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Fig. 5. — GJ 117 spectrum for the Li I 6707.8 A region plotted as a solid histogram, a. 
(top-figure) Over-plotted are PHOENIX models with 7^ ratios of 0.0 (dotted), 0.05 (dash- 
dotted), and 0.10 (long-dashed) lines. The lower panel shows x^ residuals from the difference 
of the data minus the model (see text). The interval in which x^ was calculated is indicated 
by the dotted vertical lines, b. {bottom-figure) Shows an expanded view of the lithium line 
profile with only the PHOENIX models with ^ ratios of 0.0 (dotted), 0.05 (dash-dotted) 
over-plotted, and their residuals. 
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Fig. 6. — The F statistic as a function of the 7^ ratio for Li model hne profiles as compared 
to the observations for GJ 117. F is defined as the ratio of reduced x^ for each model to 
the model with the minimum reduced x^- The F value for a 90% confidence limit (~2.80) is 
over-plotted as the dashed horizontal line. A F value of 1 is also indicated. 
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3.2. Evidence for flare activity on GJ 117 

The higher rotational velocity and deep er convection zone ra ake GJ 117 considerably 
more active than the Sun. In our earlier work (jChristian et al.ll2005l ) we presented arguments 
for the generation of Li in the atmosp h ere o f GJ 117 by spallation reactions (for example, 
see the recent work of iTatischeff et al.l (l2008l )). These arguments were based on the strong 
quiesce nt X-ray luminosity o f this source which in turn implied a high time-averaged flare 
energy (JDoyle fc Butleiill985l ). 



We have examined archival EUVE data for GJ 117 for variability. A light curve of 
these observations in the EUVE DS lOOA bandpass is shown in Figure [71 The observations 
of GJ 117 shows 3 small flares during the 540 ksec duration. The largest flare shows a 
~50% increase over the quiescent rate of 0.09 counts s~^. The variability of GJ 117 was 
also conflrmed found to be significant using the Kolmogorov-Smirnov (KS) statistical tests 
fIChristian eTaPllDDsh . 



4. Concluding Remarks 

We have obtained high resolution McDonald echelle observations of the dKl star GJ 117 



as a follow-up to our detection of ^Li in this object using VLT and UVES flChristian et al. 
20051 ). Here we report the detection of ^Li in our McDonald data at the 5% level ' ^^ 



^Li 



0.05±0.02). We have used the solar spectrum and template stars to eliminate possible blends, 
such as Ti I, in the ^Li spectral region. Additionally, bisector analysis showed no significant 
red asymmetries that would affect the lithium line profile. GJ 117 is much more active than 
th e Sun and its X - ray lu minosity is at least one order of magnitude higher. As outlined 



m [Christian et al.l (120051 ) . GJ 117 has the needed energy budget to produce the observed 
^Li in spallation reactions on the stellar surface. However, we caution that uncertainties in 
the continuum placement may cause additional errors that may decrease or increase the ^Li 
fraction for the current methods used. Future high resolution observations during a large 
fiare event on such Li-rich active dwarf stars could provide definitive proof of spallation 
reactions in these stars. 



We thank David Doss and all of the McDonald staff for their excellent support for this 
project. We thank Dr. David Lambert for suggested improvements to the manuscript. We 
also thank an anonymous referee for suggested improvements. DC thanks STFC for travel 
support. D.J. was supported by the projects No 146001 and 146007 financed by Ministry of 
Science of Serbia. 
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Fig. 7. — EUVE light curve data. The top panel shows the entire observation in 2 ksec bins 
for the Deep Survey 100 A bandpass. The lower panel shows a blow-up of the 200-400 ksec 
interval showing two flares (see text). 
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